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We present a general form of multi-dark soliton solutions of two-dimensional multi-component 
soliton systems. Multi-dark soliton solutions of the two-dimensional (2D) and one-dimensional 
(ID) multi-component Yajima-Oikawa (YO) systems, which are often called the 2D and ID multi- 
component long wave-short wave resonance interaction systems, are studied in detail. Taking the 
2D coupled YO system with two short wave and one long wave components as an example, we 
derive the general Y-dark-dark soliton solution in both the Gram type and Wronski type determi¬ 
nant forms for the 2D coupled YO system via the KP hierarchy reduction method. By imposing 
certain constraint conditions, the general Y-dark-dark soliton solution of the ID coupled YO sys¬ 
tem is further obtained. The dynamics of one dark-dark and two dark-dark solitons are analyzed in 
detail. In contrast with bright-bright soliton collisions, it is shown that dark-dark soliton collisions 
are elastic and there is no energy exchange among solitons in different components. Moreover, 
the dark-dark soliton bound states including the stationary and moving ones are discussed. For the 
stationary case, the bound states exist up to arbitrary order, whereas, for the moving case, only the 
two-soliton bound state is possible under the condition that the coefficients of nonlinear terms have 
opposite signs. 
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1. Introduction 


The two-dimensional (2D) coupled Yajima-Oikawa (YO) system, or the so-called 2D coupled 
long wave-short wave resonance interaction system: 1 


i(S^+S^)-sS+LS^=0, 

(1) 

5L 

to 

•i'T 

1 

+ 

£ 

© 

(2) 

L f = 2(o 1 |s( 1 )| 2 + a 2 |5^| 2 ). v , 

(3) 

where Oi = ±1, a 2 = ±1, S 1 ' 1 ' 1 ,e C, L,t,x,y e R, was derived as a two-component generaliza¬ 
tion of the 2D YO system (the 2D long wave-short wave resonance interaction system) by virtue 
of the reductive perturbation method. 2,3 The 2D coupled YO system can be written in the vector 

form: 

i(S, + S v )-S xv + LS=0, 

(4) 

£, = 2(||S|| 2 )„ 

(5) 

where S = (S^^S^) 7 and | | is defined as 


||S|| 2 = oi5 (1) 5 (1) * + a 2 5 (2) 5 (2) *. 

(6) 

The above 2D coupled YO system can be generalized into a multi-component system, which is 

cast into the following vector form 

i(S,+S,)-S„+LS = 0. 

(7) 

£, = 2(||S|| 2 ) V , 

(8) 


where S = <J k — ±1 for k = 1,2,• • • ,M, and || || is defined as 

M 

l|S|| 2 = ^o k S^S {k) *. (9) 

k= 1 

The one-dimensional (ID) YO system was proposed as a model equation for the interaction of 
a Langmuir wave with an ion-acoustic wave in a plasma by Yajima and Oikawa, 4 which was 
also derived from several other physical contexts. 3,5 ~ 7 The ID YO system was solved exactly 
by the inverse scattering transform method 4,8 and the (classical) Hirota’s bilinear method (which 
uses the perturbation expansion). 9,10 It admits both bright and dark soliton solutions. The 2D YO 
system for the resonant interaction between a long surface wave and a short internal wave in a 
two-layer fluid was presented and the bright and dark soltion solutions are provided by using 
the Hirota’s bilinear method. 2,3 The Painleve property for the 2D YO system was investigated 11 
and some special solutions such as positons, dromions, instantons and periodic wave solutions 
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were constructed. 11,12 For the 2D coupled case, the multi-bright soliton solutions expressed by 
the Wronskian to Eqs.(l)-(3) were provided. 1 Later, the bright N-soliton solutions in the Gram 
type determinant for the multi-component YO system were obtained. 13,14 Similar to the single 
component case, the Painleve property and dromion solutions to Eqs.(l)-(3) were discussed. 15 In a 
recent paper by Kanna, Vijayajayanthi and Lakshmanan, one and two mixed soliton solutions for 
the multi-component YO system were constructed. 16 Very recently, the rogue wave solutions for 
the single YO system in ID case were derived. 17,18 

However, to the best of our knowledge, general multi-dark soliton solutions for the multi- 
component 2D and ID YO system have not been reported yet. Moreover, general multi-dark soliton 
solutions for the multi-component 2D soliton systems have never been previously reported in the 
literature. In this paper, by using the reduction method of the KP hierarchy, we derive and prove 
the general V-dark-dark soliton solutions to Eqs.(l)-(3) and their dynamics are discussed in detail. 
Based on the KP theory, the general V-dark-dark soliton solutions expressed by either the Gram 
type or Wronski type determinant are obtained directly from the x-functions of the KP hierarchy by 
means of reductions. Similar to the ID coupled nonlinear Schrodinger (NLS) equation case, 19 it is 
very difficult to obtain multi-soliton solutions for the ID coupled YO system since some non-trivial 
constraints for parameters need to be satisfied. In this paper, we show that the general V-dark-dark 
soliton solutions for the ID coupled YO system can be obtained from the ones for the 2D coupled 
YO system by the reduction technique. Thus the constraint condition is naturally obtained. 

Kanna, Vijayajayanthi and Lakshmanan analyzed the bound states of the bright-dark solitons 16 
and the bound states of the bright-bright solitons 20 for the 2D coupled YO system. In this paper, 
we investigate the bound states of dark-dark solitons for the 2D coupled YO system. The bound 
states of dark-dark solitons for the ID coupled NLS equation with mixed focusing and defocusing 
nonlinearities was reported for the first time by Ohta, Wang and Yang. 19 The authors pointed 
out that the bound states of three or higher-dark-dark solitons do not exist. In the present paper, 
we show that the bound states of dark-dark solitons can be formed in the 2D coupled YO system 
including the stationary ones and moving ones. Lor the stationary case, the bound states of arbitrary 
order dark-dark solitons exist, whereas, for the moving case, only the bound states of two dark-dark 
solitons occur when the coefficients of nonlinear term take opposite signs. 

The rest of the paper is organized as follows. In Sect. 2, we briefly present the bilineariza¬ 
tion procedure for the 2D coupled YO system. The V-dark-dark soliton solutions with the implicit 
dispersion relation are derived through the classical Hirota’s bilinear method which uses the per¬ 
turbation expansion. In Sect. 3, the general V-dark-dark soliton solutions expressed by the Gram 
type and Wronski type determinants are obtained directly through the reduction method of the 
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KP hierarchy. Moreover, the general A-dark-dark soliton solutions for the ID coupled YO system 
are further obtained by imposing a constraint on parameters. Sect. 4 is devoted to the analysis of 
dynamics of one and two dark soltions, which suggests that the energy of solitons is completely 
transmitted through each component when two dark-dark solitons collide. In Sect. 5, the bound 
states including the stationary case and the moving case are discussed in detail. In Sect. 6, the 
general A-dark-dark soliton solutions for the ID and 2D multi-component YO system are briefly 
presented. Appendix A and B present the proofs of Lemma 2.1 and Lemma 2.4, respectively. 

2. Dark-dark soliton solutions of the two-dimensional coupled YO system 


Under the dependent variable transformation 

5 (t) = S (2) = j,L = —2(log F) xx , (10) 

Eqs. (l)-(3) can be converted into the following bilinear equations: 

[i( A +D y ) -D 2 ]GF = 0, (11) 

[i(D,+D y )-Dl]H F = 0, (12) 

{D r D x - 2 C)F ■ F + 2aiGG* 4- 2 o 2 HH* = 0, (13) 


where G and H are complex functions and F is a real function, C is an arbitrary constant and * 
denotes the complex conjugation hereafter. The Hirota’s D-opcrators are defined as 

4-A) 

2 .1 Dark-dark soliton solutions by Hirota ’s direct method 

In this subsection, we look for soliton solutions by the Hirota’s bilinear method which uses the 
perturbation expansion. 10 To this end, we expand G, H and F in terms of power series of a small 
parameter £ 

G = pi exp(i^)[1 + egi + e 2 g 2 + £ 3 g 3 + • • ■ ], (14) 

H = p 2 exp(i^ 2 )[l +e/?i +£ 2 /i 2 + e 3 /u H—], (15) 

F = 1 + e/t + £ 2 /2 + £ 3 /3 4-, (16) 

where C,j = a jx + P v y — 8/ + and p ; . a j, Py, 8y, C,jq, (j— 1,2) are real parameters. 

Substituting (14)—(16) into (11)—(13), we obtain the following constraint condition: 

C = cup 2 + G 2 p 2 , 8j = $j-a 2 j, j — 1,2. (17) 

Arranging each order of £ and solving the resultant set of linear partial differential equations re- 
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cursively, we obtain the one-soliton solution: 

G = piexp(i£i)[l +Aexp(r|)], 

H = p 2 exp(i^ 2 )[l +5exp(ri)], 
F=l+exp(r|), r\ = k x x + k y y + 03T + r\o, 


with 


(2aj k x — ky — co)i — k 2 

{ld\k x — k y — (ni)i + k 2 ’ 

B _ (2a 2 k x -k y - co)i - k 2 x 
(2a 2 fcjc — k y — co)i + k 2 ’ 

(Ok x = o lP t[2- (A+A*)] +a 2 p 2 [2- (.B + B *)], 


where k x ,k y ,co and T|o are arbitrary complex constants. 
Furthermore, we obtain the 2-soliton solution 


G — pi exp(i^i) [1 +Aiexp(pi) +A 2 exp(p 2 ) 
+Ci 2 AiA 2 exp(pi +p 2 )], 

H = p 2 exp(i^ 2 )[l +5iexp(pi) + F 2 exp(r| 2 ) 
+CnB\B 2 exp(pi +p 2 )], 

F = 1 +exp(pi) +exp(p 2 ) +Ci 2 exp(pi +p 2 ), 

with 

^ _ (2oti k x j — feyj — (Qj)i — 

{2aik x j — k y j — co 7 )i + k 2 x j 

B _ (2 CL2k x j — kyj — (0 /)i - 
( 2a2k x j — k y j — &>j ) i + 

= Pi [ 2 “ +^p] + P2 [2 - ( B j + B *j )]) 

F/ = + F/O 5 

and 


C 12 = (0i — 0 2 )(tOi —(0 2 ) +Oi P j(AiA 2 +A 2 Ai —2) 
+(52p2{ B l B 2^~ B 2 B l ~ 2 ) , 


(18) 

(19) 

( 20 ) 


( 21 ) 

( 22 ) 

(23) 
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Cl 2 — (9l + 62 ) (®1 + t 0 2 ) + °lPl (^ 1^2 +A\Al ~ 2 ) 
+o 1 pl{B l B 2 + B\B ¥ 2 -2), 


where k x j,k y j,03j and r]yo (j — 1,2) are arbitrary complex constants. 

In general, one can get the iV-dark-dark soliton solutions of the 2D coupled YO system (l)-(3): 

G = pie^ 1 

In n \ 

E ex p ILvMj + a j)+ E vw c jt , ( 24 ) 


x 


j«=o,i V=1 

// = p 2 e^ 2 


i<i</ 


x 


_ N N 

E ex p( E^(p./+( ? ./)+ E we,-/ 


/i=0,l v-i 


!</</ 


/A N 

f = E ex p E^j+ E wMji ] ’ 

p=0,i \i=l i<y'</ 


(25) 

(26) 


with 


Oipi + 0 2 p2 = C, 

£k — CL k x + — ( Pa- — oc i)t + C,ko , k= 1,2 

, aj _ a _ (2«i k x j — kyj — CQy)i — k%j 

(2a 1 k x j — k y j — toy) i + k^j 
e bj — g (2a 2 kxj~ ~ k y j — tOy)i — 

( 2a 2 k x j — k y j — toy )i + k%j 

= a iPi[2 — (Aj+A*j)\ + ct2P2[2- (Bj + #*•)], 


and 

c n 

r , = e c ;V —_ ii 

S'/ 

Cj t = (k A .y - 0/)(coy - (0,) + a lp 2(AyA; +A/A*j - 2) 
+a 2 pi(5yfi?+fi/fi}-2), 

Cjj = (. k x j + 0/) (toy + (0/) + 0i Pi (A jAj +A*jAj — 2) 
+ 02 P 2 O®/'-®/ +B*jB*i —2), 


where the notation X^ =(U represents all possible combinations /ry = 0,1 and T|y = k x jx + k y jy + 
0 )jt + riyo for j = 1 , 2,3 ... N. 
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2.2 General N-dark-dark soliton solutions in the Gram determinant form 

In this subsection, we construct an alternative form of the A'-soli ton solution based on the KP 
hierarchy reduction method. 


Lemma 2.1. The following bilinear equations in the KP hierarchy: 19 


(D X2 - T>\ - 2 aD x f)x{k + 1, Z) ■ x(M) = 0, 

(27) 

[^D x ,D x _ x -\\jX{k,l)-x{kf) 


= -x(k+l,l)x(k-l,l), 

(28) 

(D X2 - D% - 2 bD Xl )x(k,l + 1) • x(k,l) = 0, 

(29) 

(^D^Dy^-l^j x(k,l)-x(k,l) 


= -x(k,l+l)x(k,l-l), 

(30) 


where a and b are complex constants, and k and l are integers, have the Gram type determinant 
solutions 


1 <i,j<N 


i(k,l)= mij(k,l ) 
where the entries of the determinant are given by 

mij(k,l) = Cij + J§i(k,l)\\fj(k,l) dxi, 
§i(k,l) = ( pi-a) k (pi-b) l exp(Gi ), 


V/(M) = 


1 


with 


6, = 


1 


Pi-a 


-X—l + 


qt + a 

1 


qt + b 


exp(0/), 


0 , = 


1 


qt + a 


X-l + 


Pi-b 

1 

qt + b 


y~l + PiX\ + pf X 2 + 0/0, 


■y- 1 +q t x 1 - qfx 2 + 0 , 0 , 


where Cij,pi,qi , 0,o and 0,o (i,j= 1,2,--- ,N) are arbitrary complex constants. 


(31) 


The proof is given in the Appendix. 

Now we consider the reduction of the above bilinear equations in order to derive the general 
dark soliton solution. Assuming jc_i, y_i, x\ are real, X 2 , a(— iai), b(— iaf) are pure imaginary 
and q, = pf 0 /() = 0* o , c l} = cf = 5 (/ , one can get 

0 ; = 0-, m ji (k,l) = m*j(-k,-l), 
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(32) 

Hence, Eqs.(27)-(30) can be recast into 


(A 2 ~ D 2 Xl - 2i«i D xx )g-f = 0, 

(33) 

(^D Xl D x _ 1 -l)f.f=-gg*, 

(34) 

(A 2 — D 2 x — 2ia 2 D Vl )h- f = 0, 

(35) 

(^D Xl D y i — l \ f ■ f = —hh*, 

(36) 


by defining 


/ = t(0,0), g = x(l,0), h = x(0,1), 

g*=x( —1,0), h* = x(0, — 1). (37) 

By introducing the independent variable transformation 

xi =x, x 2 — -i y, 

x-i =CTipf(f-y), y-i=o 2 pl{t-y)i (38) 

i.e., 

d x = dq , dy = id^l P ] d x t ^2p2^y_i •> 


df = OlPldx-i +G2p 2 dy_ l , 

(39) 

Eqs.(33)-(36) become 


[i(A + D y — 2a\D x ) — D 2 ]g • / = 0, 

(40) 

[i(D t +D y ^2a l D x )+D 2 x ]g*-f = 0, 

(41) 

[i (A + D y - 2a 2 D x ) -D 2 x ]h-f = 0 , 

(42) 

[if A + D y - 2a 2 A) + D 2 x \h* • / = 0 , 

[AAy-2(oiPi +a 2 p2 )]/-/ 

(43) 

+2c l p 2 l gg* + 2c 2 p 2 2 hh* = 0. 

(44) 

By virtue of the following dependent variable transformation 


S^ = p 2 e^j, L=-2(log/)„, 

(45) 


with C,j — ajx + P v y — ((1/ — aj)t + C,jo for y=l,2, the bilinear equations (40)-(44) are then trans¬ 
formed into the 2D YO system (l)-(3). Hence we immediately have the following theorem for the 
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general iV-dark-dark soliton solutions of Eqs. (l)-(3). 

Theorem 2.2. The N-dark-dark soliton solutions for the 2D coupled YO system (l)-(3) are 

S l 1 ) — p ie *1 a i*+Pi:v-(Pi- a i) f +£io]£ _ 

5(2) _ p 2e i{a 2 x+$2y-{$2-a.l)t+^2o\l!_ _ 

L=-2(log/)„, 

where f,g and h are Gram determinants given by 


f = 


8 = 


h = 


Sij + ' 


1 


Pi+Pj 


-e^ 


NxN 


§ij + 


+ 


Pi - let 1 \ 1 

p*+iai J Pi + P* 


1 


with 


fi = Pj x ~ 


Pi ~ 1 CL 2 


P* + ia 2 / pi + p* 


oiPi + 


MV 


MV 


NxN 


NxN 


— ZOCi •'/ 


Pj — zaj pj — ia 2 

, f g iP? ! 


t + £>jo, 


(46) 

(47) 

(48) 


^ Pj — icL\ pj - ia 2 / 

where 0 ^. (l/ c . p/ t . ^0 (k=l,2) are rea/ constants, and PjXjo are arbitrary complex constants. 


Remark 2.3. By putting the parameters in (24)-(26) as 

= Pj+p )» = -®7 - i(p? - pf), 

_ Gipf^+p*) G 2 p^(p 7 -+pp 
^ — ioci I 2 \pj—iu 2 \ 2 

%o = ^o + ^*o + ^-2^, 

the iV-dark-dark soliton solutions (24)-(26) are equivalent to the Gram type determinant solutions 
(46)-(48). 

2.3 General N-dark-dark soliton solutions in the Wronskian form 

In this subsection, we show that the general /V-dark-dark soliton solutions for the 2D coupled 
YO system (l)-(3) can be expressed in the Wronskian from. 

Lemma 2.4. The following Wronskian satisfies the bilinear equations of the KP hierarchy (27)- 
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(30): 


with 


*(k, 1) 


9i (M) 

d Yl 9i (k,l) ■■ 

■ ar’Vw) 

92 (M) 

d Xl 9 2 (k,l) ■■ 


9 w(M) 

d*i9 a(M) ■■ 



9 i(M) = (. Pi-a) k (pi-b) l exip(Qi ) 

+ {qi-a) k {qi-b)'exp{Qj), 

n 11 9 

0/ = - X -1 H-- V- I + PiX 1 + Pi x 2 + 0/0, 

Pi-a Pi-b 

~ 1 1 2 X 

0/ =- X-1 H- -y- 1 +qiXi +cy i x 2 + 0,o, 

< 7 / - a qi — b 


where pj, c/j . 0,o and 0,o are arbitrary complex constants. 


(49) 


We provide the proof in the Appendix. 

Next, we proceed to reductions. By applying the complex conjugate conditions 

Pi = -do 0/0 = -0/0, 

with exp(0/ O ) = (nJLi.fat/gzg) exp(0' o ) and the determinant formula 21 


det lp J i Pli + qj % 

l <hj<N y 


N 


= ( A(qi,q 2 ,--- ,q N )Y[% 

i=l 


N 


X 


E I 


n 


(Piu-Pi V )(diu~di V ) 


M 

M 


V— 1 

_ l TTA Pi~9k\ m' 


=0l<ii<-<i M <N \1<ju<v<M (Pb diy'lidi/j Pi v ) 


(50) 


(51) 


where % — t B- and A is the Vandermonde determinant, the following relation can 

be derived 

N 


t(o,o) = ££ Y. 

M=01<h<-<i M <N 
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where 


Next, setting 


M 


x n ex p( 0 'v+0* v )’ 

V=1 


t(i,o ) = q-qY. E l?(v‘v)| 2 

M=0l<ii<-<i M <N 

1 N i 

<-i,o)=--§£ e i?(v<v)r 

M=01<h<-<i M <N 

m f p * + a \ 

x S(-^) exp(e '“ +e '>)’ 

T(o,i) = a-?E E l?"(v'v)| 2 

M=01<!'i<-<i M <A 

1 N i 

<0,-1) = -.?E E l?(^‘v)|’ 

M f rj* + £\ 

x S(-^) exp(e '“ +e ' J ’ 


n -^4) 

,=1 \k=\,k±i Pi Pk J 

x exp(—6*), 

G = ni- 1 )"^* + «), G = f[(-l)"(P* +4). 

i= 1 i= 1 


4(v,«v) 


n 

1</j<V<M 


(Pi,-Pi v ) 

(pi„+pi,y 


■ = #) . T(1,0) T(0,1) 

£ ’ * ’ a4 ’ 

* Gx(-1,0) Gx(0,-1) 

—-, n =-, 


(52) 


(53) 


(54) 


(55) 


(56) 


(57) 


Eqs. (27)-(30) become 


(D X2 -Dl-2aD Xi )g-f = 0, 


( 58 ) 
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(^D Xl D x _ l -l^f-f=-gg*, (59) 

(■ D X2 -D 2 Xl -2bD Xl )h-f = 0 , (60) 

(^D Xl D y _ 1 ~lJf.f=-hh*, (61) 

which are nothing but the bilinear equations (33)-(36) if a = ioci, b — ia 2 . Then, by applying the 
same transformations of independent variables (38) and dependent variables (45), the 2D coupled 
YO system (l)-(3) can be obtained. 

In summary, we obtain an alternative form of the /V-dark-dark soliton solution of Eqs. (l)-(3) 
in the following theorem: 


Theorem 2.5. The N-dark-dark soliton solutions for the 2D coupled YO system (l)-(3) are 


5 ( 1 ) — p ie 4 a i^+Pir-(Pi- a i)(+Cio]^ ; 

(62) 

S ( 2 ) = p 2 e ! '[ a 2A-+p2y-(P 2 -a5)f+C2o]^_ 

(63) 

L = -2(log/)„, 

(64) 


where f,g and h are Wronskians expressed by the form 


with 


(pl 3„(pi 


/ = 


1 

Q 


92 d*,(p2 


arv 


-JN- 1) 

(P/V o Xl cpv ••• 0*1 (P/V 



9i 

5 a, 9 i •• 

• 

1 

cp 2 

5*, 92 

■ 

ag 


5*! 9a 

■ aTV 


9i 

5v,9i •• 

■ ar^, 

i 

Cp2 

5v, 92 •• 

• 

c 2 g 

9a 

... *9- 

CD 

■ airV 


Q = \ {- p \,-p*,...-p* N ) 


A / A 

*n( n 

7=1 '£=1 Mi 


Pj + P*k \ 
P*j~PkJ 


exp (-£}), 
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and 


N 


Ci — I~[(— C) N (p*j + i‘oci), 

7=1 


N 


Ci — J^[(— l) N (p*j + i&2), 

7=1 


cp 7 = exp(y +exp(-^), 

9 y = (p ; --/ai)exp(y -(>*+zai)exp(-^*), 


9 ./ = (pj ~ /a 2 )exp(^) - + /a 2 )exp(-^), 


^j = Pi x ' 


^2P? 


P 7 
+ 


o 2 p5 


«ai pj — i a 2 


■ + ■ 


+ iPi )>’ 


^lPi 


+ ■ 


o 2 p^ 


? + ^/0> 


v Pj — uii pj — zoc 2 - 

where OCy. [)/ ( . p/ t . ^o> (&=1,2) are real constants, andPjXjO are complex constants. 


3. General Y-dark-dark soliton solutions of the one-dimensional coupled YO system 

The general TV-dark-dark soliton solutions for the ID coupled YO system can be derived from 
the one for 2D coupled YO system by further reductions. In what follows, we show the detailed 
process. 

First, it is noted that the Gram determinant solution of the bilinear equations (27)-(30) in the 
KP hierarchy can be rewritten as 


x (M) 
hj + 


_ ■-] 

•+ ia i J 


Pi- iai 


1 


Pi ~ ia 2 


P*i + ia 2 / Pi + p, 


-e^ 


= cxp ( 

,k=\ 


x 


%e-«+ 

J 1 1 rv . 1 irv« 


1 


Pj + iai 


Pi + 1«2 / 


(65) 


with 


& + 5? = 


+ 


1 


1 


■ + ■ * . 

Pi-Mi 7?-+iai 

1 1 


x-i 


+ 


y -1 


Pi - ia 2 p* + ia 2 

+ (Pi + pD*1 + {pj — p! 2 ) x 2 + ^i'0 + ^ 
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Thus, if pi satisfies the constraint condition: 

2 / 1 1 
OlPif -— + 




Pi-Ml Pj+Mi 

1 1 


+ 


Pi ~ 1CC2 Pi + !«2 


-i (ri-pf), 


i.e., 


then we have 


which implies 


sip? 


■ + ■ 


°2p2 


-i (Pi~P*), 


\pi-iai\ 2 |p;-ia 2 | 2 
(oiPi8^_! +G 2 pP y _ l )x(k,l) = -i d X2 x(k,l) , 


(oipR_! + 02 P 2 VJ/ = -i d X2 f, 

by using / = x(0,0). Moreover, differentiating with respect to xi once, we have 

Olpl/xix_i TO2P2A1V 1 = —\fx\x2■ 

Notice that Eqs.(34) and (36) can be rewritten as 

fxix- 1 f-fxi fx_ x~f 2 = -gg*, 
fx\y_\f — fx\fy_\ — f 2 = —hti* . 

From the above relations, we have 

-i(/W- fxifxx) - (aip? + a 2 p2)/ 2 

+o i p 2 l gg* + o 2 plhh* =0, 

or the bilinear form 

-i D Xl D X2 f ■ f - 2(oip? + 02 P 2 )/ 2 
+2o lP ?gg*+2CT 2 p|M*=0. 

Finally, by using the transformations of independent variables 

x\ — x, x 2 = —it, 
i.e., 


( 66 ) 

(67) 

( 68 ) 

(69) 

(70) 

(71) 

(72) 


(73) 


(74) 


(75) 


8jC] — — '0; , 


( 76 ) 
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Eqs. (33)-(36) become 


"i(A ~ 2a\D x ) — D 2 ]g • / = 0, 

(77) 

[i(D t -2a l D x )+D 2 x }g*-f = 0, 

(78) 

]i(D r -2a 2 D x )-D 2 ]h-f = 0, 

(79) 

[i(D t -2a 2 D x )+Dl]h*-f = 0, 

(80) 

D t D x -2(o l p 2 l +G 2 p 2 )}f-f 


+2o l p 2 l gg* + 2c 2 p 2 2 hh*=0. 

(81) 


By similar transformations of dependent variables (45), the above bilinear equations are converted 
into the ID coupled YO system. Thus we have the following theorem about iV-dark-dark soliton 
solutions for the ID coupled YO system. 


Theorem 3.6. The two-component generalization of one-dimensional YO system 


isl l) 


s£>+lsW = o , 


+IS (2) = 0, 

Lt = 2(G l \S^\ 2 + G 2 \sM\ 2 ) x , 


has N-dark-dark soliton solution: 


gW = p 1 ^ i [ (x i x + a i f +Cio] ^ 
— p 2e i[ a 2 x + a 2t +?2o] , 

L = -2(log/)„, 

where f,g and h are Gram determinants given by 


f = 


8 = 


h = 


8 ; y + ■ 


1 


Pi+Pj 


-e^J 


§ij + 


8/ / + 


NxN 


1 


-e^J 


Pi - ict i 

p*j + ict\ J Pi + P* 

Pi~ ia A 1 C W) 
p* + ia 2 J Pi + p* 


NxN 


NxN 


(82) 

(83) 

(84) 

(85) 

( 86 ) 
(87) 


where c^j — pjx — ipjt + ^jo, «/ 0 p/ t , ^0 (k=l,2) are real constants, pj, q v o are complex constants, 
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and these parameters satisfy the constraint conditions: 


<*iPi 


+ ■ 


G 2 p 2 2 


\Pj-iai \ 2 \p, -ia 2 


-i(Pj~P 


( 88 ) 


Remark 3.7. Compared with the two dimensional case, the parameters in the (V-dark-dark soliton 
solutions of the ID coupled YO system need to satisfy some constraint conditions. In fact, by 
rewriting the solutions (46)-(48) in the two dimensional case into the similar forms as (65), one 
can get 


& + $/ 


= (Pi + P* 


x 


<*iPi 


+ ■ 


o 2 pl 


I/?/ — ioci| 2 \pi-ia 2 \ 2 


+i (pi-p*) )y 


( OlPl <*2p| 

V|p/-iai| 2 \pi-ia 2 \ 2 


t 


+ + £*o 


(89) 


It is easy to find that the constraint conditions in the one-dimensional case are nothing but the 
zero condition for the coefficients of y in the two-dimensional case. It is interesting that the similar 
constraint conditions are obtained in finding the /V-dark-dark soliton solutions for the coupled NLS 
equation. 19 


4. Dynamics of dark-dark solitons 

4.1 Single dark-dark solitons 

To obtain a single dark-dark soliton solution in Eqs.(l)-(3), we take N 
(46)-(48). The Gram determinants read 


/=! + 


1 


Pi+P 




e = l _ 1 p l~ ia l e ^f 

Pi+P* p\ +itxi 


h = 1 


1 


p { -i«2 e ^+^ 


Pl+P* Pi +1«2 
and then the one-dark-dark soliton solution can be written as 


(1) _ P 1 ^ifocix+Pi y—(Pi —a?)f+Ciol 


S- = y e 


X 


l+is: 1 (1) -(l-^ 1 (1) )tanh 


( Si+Zj 


+ 0i 


V 


5 ( 2 ) _ P 2 c i[a ? ,.Y+p7v-(P ? ,-a;)^+Go1 

2 


1 in the formula 


(90) 

(91) 

(92) 


(93) 
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with 


x 


1 +k[ 2) - (1 — A"| 2) ) tanh 


Ai + ^ + Qi 


1 


L = ~^(pi + Pi) 2 sech 


V 2 

2 /4l +4* +®1 


(94) 

(95) 


e® 1 = 


Pi+P* 


1 

lax 


K {i) _ P l -iQCi = a\ + i(fri -QCi) 

1 pX 4“icxi a\— i(b\— ai)’ 


K (2) = P1-102 


ai +i(Z?i -a 2 ) 
/;>i+ia 2 an — i(^i — 0 C 2 )’ 
£,!+£,!= 2a ix 


2oiaipf 


+ ■ 


2o 2 aip2 


aj + (^i-ai ) 2 a] + (&i-a 2 ) 2 


+ 


2oiaipf 


■ + 


2a 2 aip^ 


j + (Z?i-ai ) 2 aj + (^]-a 2 ) 2 


-4aiZ?i y 


t + l^lOR, 


where p 1 
constant. 


ai +i^i, ai,fei,^io/f,OC/,P;,pi,^jO)(i=l,2) are real constants and ^to is a complex 



(a) 


(b) 


Fig. 1. Single dark-dark solitons (degenerate) at the fixed time t = 0 with the parameters Oi =02 = 1 ./?i = 2,pi = 
l,p 2 = 2,ai =a 2 = 0,Pi = l,p 2 = 2. 


From (93)-(95), the intensity functions of the short wave components S M ) |, | | and the long- 

2 2 

wave component L move at velocity ~ a i + ^ lP ( a ) 2 — cr^b^-a y along the x-direction. As x.y —>■ 

± 00 , |s( 1 )|^|pi|,|5( 2 )|^|p 2 |,-L^0 1 

Denoting = exp(2ip| 1 ' 1 ) and k\ 2> = exp(2i(f> j _ 1 ), the phases of the short wave components 
Si 1 ) and Si 2 ) acquire shifts in the amount of 2(|)j I ' ) and 2PJ - * but the long wave component — L 
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Fig. 2. Single dark-dark solitons (non-degenerate) at the fixed time t = 0 with the parameters c>i = 02 = l,pi = 
2,pi = l,p 2 = 2,ai = 0 ,a 2 = 2,pi = l,p 2 = 2. 


phase shifts is zero as x,y vary from — °° to +°° if 2(|)^ and 2t) ( | 2 ' 1 represent the phases of con¬ 
stants k\ 1 f and k\ 11 respectively. Without loss of generality, we can assume 2(|) ( | 1 ] , 20 1 2-! e (—7t,7t], 
(or (pj 1 \ 4>2 2 ' ) E ( — |, j]). Then the intensities of the center of the solitons (^t +4i +©i — 0) are 
|S (1) |center = |pi | cos(|)^ I) , |S (2 ) | cen ter = IP 2 I cos^f' and -L center = 2a\. For the shortwave compo¬ 
nents, the fact that the center intensities are lower than the background intensities implies these 
solitons are dark-dark solitons. 

There are two different cases corresponding to values of ai and ay. 

(i) ai = (X 2 . In this case, = k[ 2 \ i.e., <|)^ = <|)[ 2 \ this means the short wave components 
SW and s( 2 ) are proportional to each other. In this situation, the dark-dark soliton solution for 
the coupled YO system is equivalent to the dark soliton solution in the single-component YO 
system, so it is viewed as degenerate case similar to the coupled NLS equation. 19 We illustrate 
these degenerate solitons in Fig. 1. 

(ii) 0Ci ^ OC 2 . The condition k\ 1 ] ± k[~\ i.e., suggests that the components 5-^ and S^ 2) 

have different degrees of darkness at the center. In this non-degenerate single dark-dark solitons of 
the coupled YO system (l)-(3), the components S (1 J and S’ 2) are not proportional to each other. As 
is shown in Fig. 2, the intensity of the component S M ) is black, but the intensity of the component 
S {2] is gray at their centers. 

4.2 Two-dark-dark solitons 

The two-dark-dark soliton solution can be obtained by taking N — 2 in the formula (46)-(48). 
In this case, we have 


5(1) = p ie i[«i- T +Piy-(Pi- a i) f +Cio]^. 5 

fi 


( 96 ) 
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5(2) _ p 2e i[a 2 x+p 2 y-(P2-a2)f+^2o]^ ) 

fl 

(97) 

L — —2 (\ogf 2 ) xx , 

(98) 

f 2 = 1 + e ^ i +^ i + 01 - fe ^ 2+ ^ + ® 2 

(99) 

+Qi2e^ 1+ ^i + ^ 2+ ^ +01+02 , 

( 100 ) 

g 2 = 1 + tf 1 (1) e* I+ ^ +01 + K^e^ 2+ ^ +&2 


+Qi 2 K[ l) K^e^ + ^ 2+ ^ +&l+@2 , 

( 101 ) 

h 2 = 1 + k[ 2) ^ 1+ ^ +@1 +K { 2 - ) e ^ 2+ ^ +02 


+Q.\ 2 K^ e^ 1 + ^i +^ 2 +^2+® 1 +®2 ; 

( 102 ) 


Pj + P) la j 

2 _ (ai ~a2) 2 + {bi-b 2 ) 2 
(ai+a 2 ) 2 + (bi-b 2 ) 2 ’ 


& = 

Pj- iai _ 

a j + i(bj 

-at) 

j 

p*+iai “ 

cij-i{bj 

-at)’ 

k {2) = - 

Pj- i«2 _ 

aj + i(bj 

- 0 C 2 ) 

j 

p*+ia 2 

aj-Kbj 

-a 2 )’ 

Sj + %" 

— k x jx + k y jy + 0 ) j t + 2£, jQR , 


= 2a jx 

_( 2* iaj ri 2<3 2 ajp\ _ \ 

\aj + (bj — ai) 2 a 2 j + (bj- a 2 ) 2 1 / 

( 2g l fl j'Pl , 2g 2fljP2 V , 

+ (by - CC| ) 2 a 2 + (bj-a 2 ) 2 ) j0Rl 

where pj = ay+iby, ay,by,ay,[3y,py,£yo, (. 7 = 1 , 2 ) are real constants, and ^ 10,^20 are complex 
constants. 


Pi ~P2 

Z 

ci\ —a 2 + \(b\ -b 2 ) 

Pl+P* 


a 1 +a 2 + i(bi — b 2 ) 


Remark 4.8. In the case of <22 = — a\ and b 2 = bi (i.e., y ?2 = — Pp, the denominator of ^12 
becomes zero. On this critical wave number, the soliton interaction shows Y-shape. This Y-shape 
type soliton solution is called the resonant soliton solution which was found in the KP equation. 
As the two-soliton solution of the KP equation, the above two-soliton solution are classified into 
two different types of soliton interactions: 22 ' 25 
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(1) If a\d 2 < 0, 1 < ^2,i 2 - This case is called the O-type soliton interaction. In this case, two 
asymptotic soliton amplitudes 2 a\ and 2a\ (in the variable —L) can be equivalent when a 2 — 
—a \. The interaction peak (the maximum of L) is always greater than the sum of the asymptotic 
soliton amplitudes. 

(2) If ci\ci 2 > 0, 0 < Q .12 < 1. This case is called the P-type soliton interaction. In this case, two 
asymptotic soliton amplitudes 2 aj and 2a\ (in the variable — L) cannot be equivalent. The 
interaction peak (the maximum of —L) is always less than the sum of the asymptotic soliton 
amplitudes. 

Note that types of soliton interactions do not depend on the parameters b\ and /? 2 (i.e., the 
imaginary parts of p\ and p 2 ). The resonant Y-shape soliton solution is obtained by taking the 
limit b 2 —> b\ in the equal-amplitude O-type two-soliton (ci 2 = —a 1 ). The interaction coefficient 
O 12 for the two-soliton solution of the 2D coupled YO system is always non-negative although it 
can be negative for the KP two-soliton solution. 



Fig. 3. Two dark-dark solitons at the fixed time t = 0 the parameters Oi =02 = l,/?i = 1 +i,/?2 = 1.5 + i,pi = 
l,p 2 = 2,ai = l,a 2 = 2,Pi = 1,(3 2 = 2. 


The collision of two dark-dark solitons is displayed in Fig. 3. It is easy to observe that the 
two solitons pass through each other without any change of shape, darkness and velocity in both 
components after collision. Hence there is no energy transfer between the two solitons or between 
the S { 1 j and S i2] components after collision. This complete transmission of energy of dark-dark 
soliton in both components occurs not only for Gi = 02 = 1 as in Fig. 3, but also for all other Gi 
and 02 values. For the coupled YO system, this kind of phenomenon is distinctly different from 
collisions of bright-bright solitons. As reported in the paper by Kanna, Vijayajayanthi, Sakkar- 
avarthi and Lakshmanan, 13 the bright-bright solitons in the short wave components S { 11 and S' 2) 



J. Phys. Soc. Jpn. 


DRAFT 


undergo shape changing (energy redistribution) collisions while the long wave component only 
have an elastic collision. 


5. Dark-dark soliton bound states 

In this section, we investigate the soliton bound states. To obtain two dark-dark soliton bound 
states of the coupled YO system, the parameters need to satisfy and yj — which 

results in two solitons with the same velocity in both short and long wave components. 


5.1 The stationary dark-dark soliton bound states 

The stationary dark-dark soliton bound states means that the common velocity equals zero. The 
stationary solitons for the 2D coupled YO system are possible when Gi and 02 take opposite signs. 
Requiring the coefficients of t in the solution (46)-(48) of the 2D coupled YO system (01 = 1 ,02 = 
— 1) to be zero, i.e., — | p = 0. The degenerate case (CCi = a 2 ) leads to = p,, the 

dark soliton solutions for two short wave components are equivalent. The non-degenerate situation 
CCi 4 OC 2 can be further divided into two subcases ( pj = aj+ibj): 


Case (a) p\ = p^: We have bj — Kl I,** 2 and K— AT* 21 *. This case is trivial. 


42 )* 


Case (b) pf ^ p^: aj = 


2. /p?(ft.;-«l) 2 -Pi(ft.;-a2) 2 


PI-P2 


In case (b), the soliton solution is independent of the time t and | 2 — \S^ | 2 = p j — p^. Thus the 
original YO system reduces to two component linear Schrodinger equations with potential L(x,y) 
if y is viewed as the time variable. That is to say, the linear Schrodinger equation 


iS v — S xx + LS = 0. (103) 

p 2 

possess two dark soliton solutions expressed by the form (46)-(48) with the constraints p — 
\ Pi P L \2 = 0 and % = aj for k = 1 , 2 . 

Two examples of bound states are illustrated in Fig. 4 and Fig. 5, respectively. Fig. 4 shows a 
case of ^2 ^ which corresponds to an oblique bound state. Whereas, Fig. 5 displays a trivial 
case of 7 -— — which corresponds to a quasi-one-dimensional one. This kind of bound states 
can be viewed as two dark-dark soliton bound states of the linear Schrodinger equation (103) with 
potential L. 


5.2 The moving bound dark-dark soliton states 

The moving bound dark-dark soliton states require the common velocity being nonzero, i.e., 
coi 4 = 0 and Ob 4 0. Then the parameters need to satisfy the following condition: 


b\=b 2 , 
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(a) 


(b) 


(c) 


Fig. 4. The stationary dark-dark soliton bound states under the condition ^ ^ with the parameters Oi = 1, o? = 

— l,pi = y/3 + i, pi = \ + ^ 4 ^i,Pi = 1,P2 = 2,ai = l,a 2 = —2, pi = 1, (3 2 = 2. 




(c) 


Fig. 5. The stationary dark-dark soliton bound states under the condition with the parameters Oi = 1 ,02 = 

-l, pi = 1 + §i,P2 = | + §i,Pi = 1,P2 = l,OCi = 1,0C2 =2, Pi = l,p 2 = 2. 


Cl2 = 


9 

oipfa 2 - 


(flj + a ' 2 2 


+ a 2 p2« , 1 2 ( 


oip 2 (a 2 + a 2 2 ) +o 2 p 2 («? + «i 2 ) 


(104) 


where a', = b\ — ai and a 2 = fi 2 — a 2 . From the above expression, Oi and o 2 must take different 
signs. To show these moving bound dark-dark soliton states, we choose the parameters as 


c?i = — c 2 =ai=Pi=pi=&i=Z? 2 =ai=l, 

«2 = ^, ^2=2, a 2 =“, P2 = _ (1° 5 ) 

and the corresponding profiles are displayed in Figs. 6-8 at different times. 

We should point out that, in both stationary and moving bound states, the short wave compo¬ 
nents acquire non-zero phase shifts but the long wave component has no phase shift as x and y vary 
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Fig. 7. The moving bound dark-dark soli ton states with the parameters (105) at time t = 0. 



Fig. 8. The moving bound dark-dark soliton states with the parameters (105) at time t = 20. 


from —oo to +°°. This feature is the same as the general two-dark-dark solitons. Indeed, if 2(^ 1 ' 
and 2(fr 2 ' 1 represent the phases of complex constants Kj 11 and K'' ~ ] respectively, the phase shifts for 
all components are S^ ase shift = ] + 2^ 1} , S^ ase shift = 2^ 2) + 24>^ 2) and -L phase shift = 0. The 
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total phase shifts of each short wave component are equal to the sum of the individual ones of the 
two dark solitons while the phase shifts of the long wave component are always zero. 

Most recently, Sakkaravarthi and Kanna presented three bright-bright soltion bound states of 
the coupled YO system. 20 It is natural to see whether or not three- or higher-order dark-dark-soliton 
bound states exist in the coupled YO system. To ensure three- or higher-order dark-dark soliton 
bound states, at least three distinct values of pj should exist. For the stationary bound states, as two 
subcases stated in Sect. 5.1, aj can either take arbitrary positive value or is determined by bj. So 
it is not difficult to construct the stationary dark-dark bound state up to arbitrary order. However, 
for the moving dark-dark soliton bound states, from (104), all bfs values must be the same, which 
ends up at most two distinct values of aj. This observation leads to a conclusion that there is no 
three- or higher-order moving bound states. 

6. General N dark soliton solutions of the one- and two-dimensional multi-component YO 

systems 

As a matter of fact, we can extend our previous analysis to the ID and 2D multi-component 
coupled YO systems. It is known that the multi-bright soliton solutions can be derived from the 
reduction of the multi-component KP hierarchy, whereas, the multi-dark soliton solutions are ob¬ 
tained from the reduction of the single KP hierarchy but with multiple copies of shifted singular 
points. Therefore, the general dark soliton solutions for the multi-component YO systems can be 
constructed in the same spirit as the two-component YO system. The details are omitted here, and 
we present only the results for both ID and 2D multi-component YO systems. 

To seek for A-dark soliton solutions, the 2D multi-component YO system consisting of M short 
wave components and one long wave component 


i(s}* } + 4 k) ) - 4 k x +LS =0, k = 1,2, • • • ,M 

(106) 

M 

L, = 2X>|S«|J, 

k= 1 

(107) 

is transformed to the following bilinear form 


[i( A + A — 2a kD x ) — D 2 ]g£ • / = 0, 

(108) 

k= 1,2, - - - ,M 


M M 

[D,D X - 2 £ a t pj]f ■ f + 2 £ o t p \gtgl = 0, 

(109) 


k= 1 k= 1 


through the dependent variable transformations: 

= p k eAu k x+fay-{%-al)t+^ k0 ]8k, £ = 1,2 ,■■■ ,M 
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L = -2(log/)„, (110) 

where oca, Pa, Pit, Cw> are real constants. 

Similar to the procedure discussed in Sect. 2.2, taking into account the Gram type determinant 
solutions of the KP hierarchy, one can obtain iV-dark soli ton solutions as follows: 


/ = 


8k = 


5„ + — 


Pi+P, 


NxN 


&ij + 


Pi - i«A \ 1 

P*j+iUk J Pi+P*j 


+t* 

e" 3 ' a 


( 111 ) 

( 112 ) 


NxN 


with 


( M ._ ,,2 

+ E- 5£ !r'+U 

tlpj~ la k 

where pj and % jo are complex constants. 

Starting from the Wronskian solution of the KP hierarchy, iV-dark soliton solutions in the 
Wronskian form can be constructed in the same way in Sect. 2.3, which is of the following form 

(pi d*,9i • 

92 5 f| 92 • 


/= i 


aT'V. 

aTV 


9a 0-k, 9a 


s(A-l) 


J X\ 


9a 


(113) 


8k = 


with 




^i9 ( i fc) 

... a (JV_ 


1 


^9? 

... a (JV_ 

°.*i 


ag 






-(k) 

W 


... a (JV_ 
°(, 


g = 

a (-p? 

* 

,-P2N- 

) — Pa) 



A / 

*n 

' A 

n 

Pj + P*k\ 

exp(- 


7=1 ' 

( A= 1 ,k±j 

Pj-PkJ 



(114) 


A 


Ck — I~[( — l) N (p*j + l a it)) 

7=1 



J. Phys. Soc. Jpn. 


DRAFT 


and 


(p; = exp(y + exp(— 


<Pj= (P j - m) exp (y - (p* + ice*) exp (-£*), 


%i = PF 


M 


Gkpj 


Pj-m 


+ i p 2 i )y 


M r- n 2 

+ E^7T'H;o, 

n p.i - ia t 


where pj and q 7 o are complex constants. 

By the similar procedure discussed in Sect. 3, TV-dark soliton solutions for ID multi-component 
YO system are provided with the same form as the one for two-dimensional case. In other words, 
TV-dark soliton solutions for ID and 2D integrable systems can be deduced simultaneously without 
reformulating the problem. To be more specific, the following bilinear form 


[i(A - 2a k D x ) - Dl\g k ■ f = 0, k = 1,2, • • • ,M 

M M 

[P,D, - 2 £ GtfDf ■ f +2 £ atfilsal = 0, 

k= 1 k= 1 

is converted from the ID multi-component YO system 

iS, ( * } - S { $ +LS® = 0, k = 1,2, • • • ,M 

M 

U = 2£o t |S<‘)|J. 

k= 1 

through dependent variable transformations 


S (fc) = p bjlukx+a-lt+fa}^ 

d 2 

L=_2 d? l0g/ ’ 


where c^, p k, are real constants. 

It is shown that, by imposing the constraint conditions 


M 


°*P* 


\Pj~ ia k\‘ 


= “i (Pj 


■Pi) 


j=l,2-,N, 


(115) 

(116) 

(117) 

(118) 

(119) 

( 120 ) 


( 121 ) 


the terms associated with D y in (108)-(109) are dropped out, the bilinear equations for the 2D multi- 
component YO system are reduced to (115)-(116) for the ID case. Therefore, the TV-dark soliton 
solution for the ID multi-component coupled YO system shares the same Gram determinant form 
(111)-(112) or Wronskian form (113)-(114) except the constraint conditions (121). 
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7. Conclusions 

We have constructed the general multi-dark soliton solutions in both the ID and the 2D multi- 
component coupled YO systems and analyzed their dynamical behaviors. General multi-dark soli¬ 
ton solutions for 2D multi-component soliton systems have never been reported in literature. 

By using the classical Hirota bilinear method, we have presented the Y-dark-dark soliton solu¬ 
tions with the implicit dispersion relation in the 2D coupled YO system containing two short wave 
component and one long wave component. By virtue of the reduction method of the KP hierarchy, 
Y-dark-dark soliton solutions expressed by Gram type and Wronski type determinants are derived 
and proved. The process of obtaining Y-dark-dark soliton solutions elucidates the connections of 
the YO system with other integrable systems in the KP hierarchy, which will be helpful for the 
further study of these systems. By further reduction, we also provide the general Y-dark-dark soli¬ 
ton solutions for the ID coupled YO system in the same form as the one for the 2D coupled YO 
system except some constraint conditions. The similar form of general Y-dark soliton solutions in 
the ID and 2D multi-component YO systems are constructed by simply inserting more copies of 
the shifts of singular points. 

We have further investigated the dynamical behaviors of one and two dark-dark solitons in the 
2D coupled YO system with two short wave components. In contrast with bright-bright soliton 
solutions, it is shown that dark-dark soliton collisions are elastic and there is no energy exchange 
in two components of each soliton. 

Moreover, the dark-dark soliton bound states including the stationary and moving ones are 
discussed. For the stationary case, the bound states exist up to arbitrary order, whereas, for the 
moving case, only two-soliton bound state is possible under the condition that the coefficients of 
nonlinear terms have opposite signs. 
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Appendix A 

In this appendix, we present the proof of Lemma 2.1 in Sect. 2.2. 19 Consider functions (]),■ and 
\\fj which satisfy the following differential and difference rules: 

= tyi(k- 1,0, 

tyi(k+ 1,0 = (3 Xl -a)<|>i(M), (A-l) 

3* 2 V*(M) = Vi(M)> 

9^_ 1 V,-(fc,0 = -Vi(fc+1,0> 

= -(d Xl +a)\|//(fc,0. 


Define 

mij(k,l ) = c,y + y <]);(£, /)\|/y(A:,/) d.vi, 

and the N x N matirx m(k,l) = (nijj(kj)) i<;j<jv. Then one can easily verify that the matrix ele¬ 
ments nijj(k,l) satisfy 

d Xl niij(k,l) = §i(k,l)\\fj(k,l), 
d X2 nijj(k,l ) 

— (d Xl tyi(k,l))\\fj(k,l) <f) 7 (/r,/) (d Xl \\fj(k ,/)), 
dx-i*nij(k,l) = —tyi(k—l,l)\\fj(k+l,l), (A-2) 

m,- ; -(k + 1,/) = mjj(kj) +§i(k,l)\\fj(k+l,l ). 


The functions (J)/, \j/,- and the matrix elements nijj(k,l ) in Lemma 2.1 satisfy these relations. 

Then with the help of (A-l) and (A-2), one can check that the derivatives and shifts of the x 
function are expressed by the bordered determinants as follows: 10,19 


d X \ ^{k, /) — 


m(M) <*>(M) 
-¥(£, 0 o 


d 2 Xl x(k,l) = 


m(M) d x MkJ) 

- V L( k,l ) 0 


+ 


m (k,l) <Z>(k,l) 

d Xl } V(k,l) 0 



d X2 x{kJ) = 




m (k,l) o 

»*iMM) 

-MM) 

0 

m (k,l) 

MM) 

-3* MM) 

0 

m (k,l) 

<P(k-l,l) 

^(k+l ,/) 

0 


-1 MM) 

m{k,l) <P(k-l,l) <t>(k,l) 
^(k+1,1) 0 -1 

-¥(*,/) -1 0 


x(k+l,l) = 

x(k-l,l) = 


m(k,l) <P(k,/) 

-V(k+l,l) 1 

m (k,l) <P(k-l,l) 

MM) 1 


(3*, + a)x(k +1, /) — 


m (M) d Xl <t>(k,l) 

- x ¥(k +1,1) a 


(d Xl + a) 2 x(k + 1,/) 
m (k,l) 

+ ~^(k+l,l) 

-'F(M) 


m (M) d 2 Xl <t>(k,l) 

- x ¥{k+1,1) a 2 

d Xl ®i(k,l) M(M) 

a 1 , 

0 0 


(d X2 + cr)x(k +1,1) 
m (k,l) 

- -^{k+1,1) 

-^(k, l) 


m (M) ^(M) 

-^(k+1,1) 

MMM) MM) 

a 1 

0 0 
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where the bordered determinants are defined as 



m (k,l) 

*(M) 



- 

-nk,l) 

0 




mn(k,l) 

mn{k,l) 

m\ N {k,l) 

4>1 (*,/) 


m 2 i(k,l) 

m 22 {k,l) 

m 2N (k,l) 

Hk,l) 


m N \ (k,l) 

m N2 (k,l) 

m NN (k,l) 

§N{k,l) 


-Vi (k,l) 

-V2 {k,l) ■■ 

■ -VN(k,l) 

0 


and ®(k,l) = (<|>i(M)><I> 2 (M),-" ,<|>a(M)) T and'?(£,/) = (\|/i (k,l), \|/ 2 (k ,/),•• • ,\\r N (k,l)). 
By using the above relations, one can verify 


(D X2 - D 2 Xl - 2aD X] )x(k + 1, /) • x{k, l ) 

m (M) ^(M) <*>(M) 

= —2 «-¥(*+1,Z) a 1 

-^(M) 0 0 

x m(fc,/) 

+2 m(M) d x MW) 

-^(k+1,1) 0 

m(k,l) <P(k,l) 

-'F(M) 0 

_ 2 m (k,l) &(k,l) 

-^(k+1,1) 0 

m (k,l) d Xl &(k,l) 

x 

-'F(M) 0 


(A-3) 


(jD Xl D x _ x - l) x(k,l)-x(k,l) +x(k+l,l)x(k- 1,/) 

m(M) <!>(£-1,Z) <D(ik,Z) 

^(k+lj) 0 -1 

~'¥(k,l) -1 0 

x m(k,l) 
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m(M) <fr(M) 

-¥(*,/) 0 

m (k,l) <Z>(k-l,l) 

x 

^{k+1,1) 0 

m(M) <J>(M) 

- x ¥{k + 1,/) 1 

m (k,l) ®(k-l,l) 

x v v . (A-4) 

¥(*,/) 1 

Both (A-3) and (A-4) are identically zero because of the Jacobi identities 10 and hence T(kJ).x(k + 
1, /) and x(k —1,1) satisfy the bilinear equations (27) and (28). In a similar way, one can prove the 
other two bilinear identities (29) and (30). 


Appendix B 

Here we present the proof of Lemma 2.4 in Sect. 2.3. 

Consider the x-function 

x (M) 

9i(M) djqtpi(M) ••• ’Vt(M) 
(P2(^,0 3jq(p2(^,0 •" air l) Cp 2 (^/) 

<?N(,k,i) a Xl cp N {k,i) ■■■ a x f ^cpa^m) 

where functions cp,(A:, /) satisfy the following linear dispersion relations: 

a,-, 9/ n) (M) = ( P/ ,l) (k + 1, /) + acp ; -' !) (fc, /), 

a X2 c Pi (n) (c/) = aj cp ; (n) (c/) 

— ^ (fc + 2,/) + 2acp ; - 1,/)+<7 _ cp ; - \k,l), 

d x _^ n) (kj) = ^ ) (k-\,l), 

3*, <PS M) (A:,/) = q>f° (M + 1) + (A, 0, 

a X2 cp 7 w (^/) = a“ 1 cp ; ( ' ,) (fc,/) 


(B 1) 

(B-2) 


(B-3) 

(B-4) 

(B-5) 
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= <pS"’ (M + 2) + 2fcip)" ) (*,/ + 1) + b 2 <f l l " > ( k,l ), 

The functions cp j(k,l) in Lemma 2.4 satisfy these relations. 

Let us introduce a simplified notation, 

\nk,l,n + l kj i,...,n+N- 1 k ,i\ 

^<piCU) (*,/) ••• 

<92(L0 a" +1 cp2(*,/) ••• d n + N -^ 2 (k,i) 

^VOM) ••• 

One can rewrite the above T-function as 

x(k,l) = \Qkj,lk,i, ■ ■ ■ ,N — 2kj,N — lk,i\ 

= |0 k,h 4,Z> ■■■jN — 2k ; i, (N - 2 k+ ij) +a(N - 2 k j) | 

= |0jt,z, U,/,.. .,N-2 k j,N-2 k+ ij\ 

= \Ok,hOk+i,h U+t ,ii ■ ■ ■ ,N— 3 k +ij,N — 2fc+i,/1 

= |0jt,/,0jfc + i j /,0jt + 2 I /, • • •,Ofc+JV—2,/,Ofc-i-jv— 1,/1 • 

For simplicity, we omit subscripts k j. Thus the above x-function is written as 

x(M) = |0jt,/, 4,/,... ,N — 2kj,N~ l k j\ 

— |0jfc,/)0fc+i,Z)0^+2,/, • • •, Ofc+jv—2,; 5 Ofc+jv— l,/ 1 

= |0,1,... ,N — 2,N — 1|. 

The differential formulas for x(k, l ) are derived as follows: 

a^T(fc,/) = a^|o,i,...,iv-2,iv-i| 

= \0,1,...,N-2,N\+Nax(k,l), 

3?,x(/c, /) = 3 X1 10,1,... ,iV - 2,N\ +Nad Xl x{k , /), 
= |0,l,...,iV-3,iV-l,iV| 


(B-6) 

(B-7) 


(B-8) 


(B-9) 
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+ \0,l,...,N — 3,N — 2,N+l\ 

+Na(d Xl — Na)x(k,l) + Nad Xl x(k,l) 

= \0,1,...,N-3,N-1,N\ 

+ \0,l,...,N-3,N-2,N+l\ 

+2 Nad Xl x{kJ) — N 2 a 2 x(kJ), (BIO) 

a X2 x(M) = a* 2 |o,i,...,iv- 2 ,jv-i| 

= — \0,l,...,N — 3,N—l,N\ 

+ \0,l,...,N-3,N-2,N+l\ 

+2 ad Xl z(k,l)-Na 2 z(k,l), (B-ll) 

d x _ l x(k,l) = d x _ l \0,l,...,N-2,N-l\ 

— | — 1,1,...,AA — 2,2V — 1|, (B-12) 

^x- 1 9xi 0 

= dx-i |0,1,... ,N — 2,N\ +Nad x _ 1 x(k,l) 

— | — 1,1,...,AA — 2,2V| H- — 2,AA — 1| 

+Nad Xl x(kJ) 

= \^X,l,...,N-2,N\ + x(k,l)+Nad x _ l x(k,l). 

(B-13) 


Consider the following determinant identities: 


1 

N — 2 

N- 1 

N 

N + 1 


0 


0 


N—l 

N 

N+l 

0 


N — 2 


and 


\ 

-i ; 

l 

N-2 

N 


0 


N- 1 

-l | 

0 


N 

0 


N-2 

N- 1 


Applying the Laplace expansion to the left-hand side of these identities, we obtain the Pliicker 
relations 


11,. - - ,N—2,N,N+l\ x |0,1,...,A-2,A-1| 

— 1 1, ... ,N—2,N—l,N+l\ x |0,1, ... ,N—2,N\ 

+ 1 1, ... ,N—2,N—1,N\ x |0,1,...,A-2,A+1| =0, 
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+ \l,2,...,N-*2,N-l,N\ x |-l,0,l,...,iV-3,JV-2| =0. 
By using the T-functions, these determinant identities are rewritten as 

2 ( — ^ 2 X (^ +1)0 + ^ X| x(A: +1,/) 
-2(N-l)ad Xl x{k+ 1,0 + iV(iV- l)a 2 x(fc+ 1,/)) 
xx(k,l) 

— (d X [X(fc+ 1,/) — Nax{k-\- 1,0) 
x(3 Xl x(M) - Nax(k,l )) 

+x(fc +1, /) 

x 1 (a X2 x(fc, /)+a 2 x(fc, /) - 2(iv+ i)aa xi x(/:, /) 
+iV(iV + l)a 2 x(fc,/)) =0, 

(a x _ i a xi x(fc,/) —x(k.l) — Nad x _ l x(kJ)) x x(M) 
-a x _jX(fc,/) X (3 Xl x(fc,Z) — Nax(k,l)) 

+x(& + 1,0 x x(fc —1,0 = 0, 


which are nothing but bilinear equations (27) and (28), respectively. Eqs.(29) and (30) can be 
proved in a similar way. 
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